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ABSTRACT 

The impact of Type la supernova ejecta on a helium-star companion is investigated via high- 
resolution, two-dimensional hydrodynamic simulations. For a range of helium-star models and initial 
binary separations it is found that the mass unbound in the interaction, (JMub, is related to the initial 
binary separation, a, by a power law of the form 5M^\j cx a™. This power-law index is found to 
vary from —3.1 to —4.0, depending on the mass of the helium star. The small range of this index 
brackets values found previously for hydrogen-rich companions, suggesting that the dependence of the 
unbound mass on orbital separation is not strongly sensitive to the nature of the binary companion. 
The kick velocity is also related to the initial binary separation by a power law with an index in a 
range from —2.7 to —3.3, but the power-law index differs from those found in previous studies for 
hydrogen-rich companions. The space motion of the companion after the supernova is dominated 
by its orbital velocity in the pre-supernova binary system. The level of Ni/Fe contamination of the 
companion resulting from the passage of the supernova ejecta is difficult to estimate, but an upper 
limit on the mass of bound nickel is found to be ~ 5 x 10"'' M©. 
Subject headings: binaries: close — supernovae: general — methods: numerical 



1. INTRODUCTION 

Type la supernovae (SNe la) are among the most 
catastrophic and energetic events in the Universe. Be- 
cause they are so luminous, and because their light- 
curve shapes and absolute magnitudes are correlated, 
SNe la can be used as "standardizable candles" in mea- 
suring the distances to remote galaxies, allowing con- 
straints to be placed on key cosmol ogical parameters 
(jBranch fc Tammannlll992l : lRiesdll996f ). SNe la also play 
a major role in galactic chemical evolution via their en- 
ergy and metal input to the interstellar medium. Thus, 
the nature of their progenitor systems and the physi- 
cal origin of variations in their properties are fundamen- 
tal problems of g reat interest (Domfnguez et al. 200l|; 
IHowell et al. l [2007l) . 

Most work on progenitor models for SNe la has fo- 
cused on two general classes of systems. The single- 
degenerate scenario involves a CO white dwarf (WD) ac- 
creting matter from a non-compact stellar binary com- 
panion, event ually becoming unstable to explosiv e nu- 
clear burning (|Whelan fc Ibenlll973t [Nomotolll982h . On 
the other hand, the double-degenerate scenario involves 
the merger of two CO WDs whose orbital decay results 
from the loss o f angular momentum due to gravitatio nal 
wave emission (llben fc Tutukovl[T98llWebbinkl[T98l . 

Single-degenerate models require that the rate of mass 
accretion be such that the WD avoid a nova explosion. 
For a high accumulation efficiency, stable burning on 
the WD is required, but this burning phase occurs in 
a fairly narrow range above 10~^ Mq yr~^. For rates 
greater than this range, a wind from t he WD devel- 
ops that limits the accretion efficiencv ([Nomotol 119821 : 
IHachisu et all 119961: llvanova fc TaamI |200J), making it 
difficult to explain the observed low-redshift SN la rate 



of 3 X 10"'^ Mpc ^ yr"i (iMannuccil |2005[) using single- 
degenerate progenitors alone. Although the predicted 
and observed numbers of double-degenerate syste ms are 
suffic i ent to explain the observed ra te ([Neleman s et al.l 
[200lt INapiwotzki eTal] [200ll [20021 ). it is more likely 
that WD- WD mergers lead to the production of an 
ONeMg WD fo llowed by accretion-iiiduced collapse to 
a neutron star (iNomoto fc IbenI [ 1985': 'i vanova fc TaamI 
120041: IDessart et al.ll2006l: IWickFa masing he et al.ll2009r 
Moreover, the wide range of total WD masses suggests 
that double-degenerate models should exhibit much more 
heterogeneity than is observed in SN la light curves 
(IGoldhaber et al.l[200irKnoD et al.l[2003l) . 

iMaod ( 2008[) arialvzed the observed SN la rate for as- 
sumed initial mass functions (IMF) and concluded that 
almost all intermediate-mass close binary systems in the 
range 3 — 8 Af© should evolve to the SN la stage. This 
result allows for a wide range of possibilities for SN la 
progenitor systems. Furthermore, observations of the 
SN la rate as a function of redshift suggest the need for 
a two-component model for the delay time distribution 
(DTD). (The delay time is the interval between a star's 
arrival on the zero- age main sequence and i ts des truc- 
tion in an SN la.) IScannapieco fc BildsteiJ (j2005D and 
iMannucci et al.l (|2006l) found that the observations can 
be fit with a short-delay-time population having delays of 
^ 10^ yr and a long-delay-time population having delays 
of 3 - 4 Gyr. 

To provide an explanation for the two populations 
suggested by the observed SN la rate, several pre- 
supernova progenitor models have been investigated. 
The long-delay-time population can be understood 
in terms of progenitor systems characterized by a 
main-sequence-like companion in the MS- WD channel 
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(|Hachisu et all l2008bD and/or by a red gia nt in th e 
RG-WD channel (iHachisu et al.l[l996l Il999allbl . [2008al) . 

In contrast, the short-delay-time population may con- 
sist of systems characterized by a massive MS star 
in the MS- WD channel or by a He star i n the He- 
WD channel (IWaldman et al.l[200llWang et al. 2009a^b|; 
IWang fc Hanll2010l:IMeng fc Yandl2009D . 

Numerically. [Marietta et al.l ( 200C1[) explored the in- 
fluence of the supernova explosion on the companion 
star in the single-degenerate channel with hydrogen- 
rich stars consisting of red giants, subgiants, and main- 
sequence stars using two-dimensional Eulerian hydrody- 
namics simulations. They found that significant quanti- 
ties of hydrogen would be unbound from the companion 
star in each case (15% of the envelope for main-sequence 
and subgiant cases, and 98% of the red giant envelope), 
in conflict with observational upper limits on the amount 
of hydrogen inferred from SN la spectra (jMattila et al.l 
120051 ILeonar"dl[2007h . 

A m ore recent hydrodynamical study bv lPakmor et all 
(120081) reexamine d the main-sequence simulation 
of [Marietta et al.l (|2000D using a three-dimensional 
smoothed particle hydrodynamics (SPH) simulation. In 
contrast to Marietta et al., Pakmor et al. adopted the 
structure for the comp anion star based on the binary 
evolutionary models of llvanova fc Taamj ()2004D , which 
yielded more compact main-sequence-like companions. 
As a consequence, Pakmor et al. found a tenfold reduc- 
tion in the amount of unbound mass compared with 
Marietta et al., bringing the prediction of the amount 
of unbound hydrogen-rich material into agreement with 
the observational upper limits. 

In a complementary analytical study, iKasenI ()2010[ ) in- 
vestigated the radiation emitted by the collision of SN la 
ejecta with a red giant, finding that the light curve should 
depend on the viewing angle. This result was attributed 
to the fact that the gas is more transparent in the re- 
gion shadowed by the companion star. This suggests the 
possibility that the secondary star may be detectable in 
future observational studies. 

These previous simulations of the effect of a super- 
nova impact on a companion star have been carried out 
only for models applicable t o the long-de l ay-time p opu- 
lation of SN la. In contrast. IWang et all ()2009allbl ) sug- 
gest a progenitor binary model based upon a helium-star 
channel for the short-delay-time population. Using a bi- 
nary evolution model, for this channel they find a SN la 
birthrate ~ 3 x 10~^ yr^^ and a corr esponding delay tirn e 
of - 4.5 X 10^ yr to - 1.4 x 10^ yr ()Wang fc Hanll20ll . 
The latter delay time is consistent with that estimated 
from observations for the short-dclay-tim e population b y 
[Scannapicc o & Bildstcn (2005), Mannuc ci et al.l (|2006f ). 
and lAubourg et al.l (|2008l ). 

In this paper we report the results of Eulerian hy- 
drodynamics simulations of the impact of SN la ejecta 
on companion stars for the single-degenerate helium-star 
channel. In the next section, the assumptions underlying 
our study, the construction of the initial model, and the 
numerical method are described. Our numerical results 
for a range of helium-star models and orbital separations 
are reported in § 3 and their implications are discussed 
in § 4. In the final section, we summarize our results and 
make some concluding remarks. 
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Fig. 1. — Density profiles of the helium-star models at the onset 
of the supernova explosion, as obtained using the onc-dimensional 
stellar evolution code. 

2. NUMERICAL METHODS AND MODELS 
2.1. Numerical codes 

For our hy drodynamical simula t ions we used FLASH 
version 3 (|Frvxell et al.l I2OOOI: IDubev et al.l I200I . 
FLASH is a parallel, multi-dimensional hydrodynam- 
ics code based on block-structured adaptive mesh re- 
finement (AMR). To solve the Euler equations on the 
AMR grid, we used the p iecewise parabolic method 
(jColella fc Woodwardl[T98^ with modifications t o han- 
dle nonideal equations of state (jColella fc GlazI [1985) . 
The equation of state used is interpolated from a precom- 
puted table of the Helmholtz free energy. It includes con- 
tributions from radiation, completely ionized nucle i, and 
degen erate electrons and positrons (Timmcs & Swestvl 
I2OOOD for an optically thick mixture of gas and radiation 
in local thermodynamic equilibrium. 

The helium-star models used in our simulations were 
generated using a one-dime nsional stellar evolution code 
(|Eggletonl[T97l[ 119721 [l973l ). We simulated four helium- 
star models with initial masses equal to 1.25, 1.35, 1.4, 
and 1.8 Mq. To evolve the helium-star models to the on- 
set of the supernova explosion, an artificial constant mass 
loss rate was adopted such that the evolution time and 
final helium star masses were c onsistent with the de tailed 
binary evolutionary models of IWang et al.l (|2009b( ). The 
resulting models are summarized in Table [1] The den- 
sity profiles of the helium-star models at the onset of the 
supernova explosion are illustrated in Figure [1] It can be 
seen that the more massive models are characterized by 
larger radii and less compact cores. 

2.2. Initial setup 

Since the speed of the ejecta in a SN la (< 10* km s^^) 
is much higher than the orbital speed of the helium star 
(< 10"^ km s""'^) in a binary system, we ignore the or- 
bital motion in the first approximation and consider a 
2D axisymmetric geometry. The simulation domain is 
described using cylindrical coordinates {r,z), with the 
z-axis defined as the direction along the line connecting 
the centers of the white dwarf and the helium star. We 
consider a simulation domain with a size equal to fifteen 
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TABLE 1 
Helium-star models. 





Mass 


Radius 


Evolution time 




Cub^ 


nikick 


'-'kick 


Model 


(Mq) 


(10^° cm) 


(yr) 










He-WDa 


0.697 


0.63 


9.2 X 10" 


-4.01 


1.42 


-3.28 


4150 


He-WDb 


0.803 


1.10 


5.2 X lO'"^ 


-3.13 


0.70 


-2.90 


2413 


He-WDc 


1.007 


1.35 


2.2 X lO'"^ 


-3.48 


1.17 


-3.18 


2703 


Hc-WDd 


1.206 


1.63 


6.1 X lO"' 


-3.51 


1.35 


-2.71 


1729 



"The entries for the power-law index, mub, and power-law constant, Cub, refer to 
the power-law fit to unbound mass versus orbital separation described in eq. [T] 

''The entries for the power- law index, miji^k, and power-law constant, Ckick, refer 
to the power-law fit to kick velocity versus orbital separation described in eq. [2] 



times the radius of the hehum star (i?He) in the radial 
direction and 30 i?Hc in the axial direction. For con- 
venience, the helium star is located at the origin of the 
coordinate grid. 

To simplify the problem, the composition of the one- 
dimensional helium-star model was taken to be a uni- 
form distribution of 98% helium and 2% carbon by mass 
when used in FLASH. This simplification leads to an 
error of < 2% in the composition and < 4% in the ra- 
dius for the lowest-mass helium-star model. To initialize 
the two-dimensional FLASH simulations, we first inter- 
polated the one-dimensional model onto the FLASH grid 
using up to twelve levels of refinement based on the mag- 
nitudes of the second derivatives of gas density and pres- 
sure. With each block containing 8 x 16 zones, the equiv- 
alent uniform-grid resolution is thus 16, 384 x 32, 768. 
For model He-WDc, the minimum zone spacing at this 
level of refinement corresponds to 1.22 x 10^ cm. Within 
the helium star, a minimum of nine levels of refinement 
was used (corresponding to a maximum zone spacing of 
9.76 X 10^ cm for model He-WDc). Because the surface of 
the helium star is characterized by a very steep density 
gradient, we established a sharp cutoff radius to avoid 
pressure errors at the surface. The cutoff radius was 
chosen such that each density drop of one order of mag- 
nitude is resolved by at least three zones in the surface 
region (we have more than 40 zones per order of magni- 
tude in the core region). Outside the cutoff radius, the 
density was set to an ambient value of 2.25 x 10~^ g cm~^, 
and the pressure was set to the value of the helium-star 
model pressure profile at the cutoff radius. The simula- 
tions, therefore, employed the same ambient density but 
somewhat different ambient pressures. The composition 
of the ambient gas in each case was taken to be pure 
hydrogen. 

The helium-star model was relaxed on the Eulerian 
grid by artificially damping the momentum for a period 
of time greater than about thirty times the average dy- 
namical time scale. During this time the damping fac- 
tor was smoothly increased from 0.7 to 0.99, ensuring 
that the Mach number in the helium-star interior was 
always smaller than 0.01. Once this process was com- 
plete, the damping was removed and the gas velocity 
was reset to zero. The supernova explosion was then in- 
troduced to the grid. During the subsequent evolution, 
we allowed second-derivative refinement up to seven lev- 
els everywhere except in two regions: within the helium 
star at any time, and in a region surrounding the su- 
pernova explosion for the first 150 seconds, we required 



nine levels of refinement (equivalent to a resolution of 
2048 X 4096 for a uniform grid with a minimum zone 
spacing of 9.77 x 10^ cm in model He-WDc). The extra 
refinement for the supernova region reduces the influence 
of grid artifacts on the developing explosion. The explo- 
sion itself was introduced by creating a spherical region of 
high-density and high-temperature gas with radius equal 
to twenty times the minimum zone spacing at nine lev- 
els of refinement. Each run used this "7/9" refinement 
pattern except for one model that was also run at 6/8, 
8/10, 9/11, and 10/12 to study convergence (see § I3.3|) . 

The Type la supernova model used is the W7 model 
described bv lNomoto et al.l ()1984f ). which corresponds to 
a carbon deflagration in an accreting CO WD. The initial 
carbon deflagration is assumed to develop in the central 
regions of the white dwarf and is then followed by a det- 
onation. This one-dimensional model provides a good 
fit to observed light curves and can be approximated by 
a white dwarf of mass Afwd = 1-378 Mq, total explo- 
sion energy E^n = 1.233 x 10^^ erg, and average speed 
Vsn — 8.527 x 10^ km s"^. This mass is uniformly dis- 
tributed within the spherical perturbation used to start 
the calculation. The velocity inside the perturbation is 
taken to be radially outward and uniform in magnitude; 
the internal energy is set using Vsn and Egn- Wc assume 
the ejecta to be entirely ^^Ni and use this fiuid compo- 
nent as a tracer for the ejecta. The mass of the supernova 
ejecta creates a potential perturbation that changes the 
equilibrium state of the helium star. However, in all of 
our models the radius of the helium star is smaller than 
the Roche lobe radius, and the timescale on which the 
equilibrium is upset by the supernova shock, a/vsn (a 
is the initial binary separation) , is shorter than the aver- 
age dynamical timescale in the helium star. Thus we ex- 
pect that the altered pre-supernova potential should not 
have a significant effect on the structure of our models. 
However, because the orbital motion is ignored in this 
study, the gravitational force from the supernova mate- 
rial will attract the helium star and cause a small velocity 
perturbation toward the supernova material. This veloc- 
ity perturbation can be comparable to the kick velocity 
of the helium star after supernova impact but is much 
smaller than the ejecta speed. 

3. RESULTS 

In this section, we describe the numerical results for 
the standard case, model He-WDc, and explore the de- 
pendence of the system's evolution on the mass of the 
helium-star companion and the initial binary orbital sep- 
aration a. To determine the sensitivity of the results to 
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the numerical resolution, we also describe a eonvergenee 
test. 

3.1. Qualitative description of evolution 

Immediately after the onset of the supernova explo- 
sion, a double shock structure is formed as the ejeeta 
interact with the surrounding medium (label A in Fig- 
ure [21 the figure illustrates the evolution of a simulation 
with initial binary separation of 4 x 10^° cm, ^ 3 i?Hc, 
and 9/11 levels of refinement). A forward shock ex- 
pands outward into the ambient medium, and a reverse 
shock propagates inward (in Lagrangian coordinates). 
The two shocks are separated by a contact discontinuity 
(|Dwarkadas fc Chevalieii[T998D . During this free expan- 
sion phase, the swept-up ambient medium has very little 
effect, except for the development of an instability at the 
contact discontinuity. This instability is unimportant in 
affecting the helium star, as it is seeded by a numeri- 
cal grid effect and is Rayleigh- Taylor unstable (label A) 
(|Dwarkadaj2000[ ). The expansion of the ejecta leaves be- 
hind an extremely low-density region at the center of the 
supernova (label a) that is prone to numerical grid effects 
which become visible due to the color scale; these are 
smoothed out later by the reflected shock (label C). The 
SN ejecta reach the companion at i ~ a/vsn ~ 50 sec- 
onds, at which time a bow shock forms at the leading 
surface of the companion, making an angle of 40° with 
respect to the z-axis. This can be seen in the upper- left 
panel of Figure [2] (label B^Q. As the bow shock prop- 
agates further, the shearing of gas in conjunction with 
the action of gravity due to the helium star causes dis- 
tortions in the bow shock structure. After the impact 
of the ejecta on the helium star, ejecta material begins 
to reflect and refill the central supernova region (label 
C in Figure [2]). As a result of the mixing beginning at 
^100 seconds associated with the Kelvin- Helmholtz and 
Rayleigh- Taylor instabilities in the envelope of the he- 
lium star (label D in Figure [5]), the bow shock breaks up 
and is divided into a curved shock and a straight shock 
that makes an angle of ~ 40° with the z-axis (label B 
in the upper- right panel of Figure [2]). The instabilities 
continue to develop between ~ 100 s and ^ 600 s as the 
mixing region moves away from the helium-star envelope 
(lower-left panel of Figure^]). At around 100 s, the shock 
propagates to the center of the helium star, leading to the 
pulsation of the helium star. Smooth shocks (label E in 
Figured]) are continually generated by the oscillation of 
the helium star after ^ 600 s. Radial shocks (label F in 
Figure [5]), which result from the interaction of shocks on 
the rear side of the helium star, sweep around the he- 
lium star (lower-right panel of Figure [2]) . As the helium 
star evolves further, its density profile becomes smoother, 
approaching a new equilibrium state. Q ualitatively, our 
results are similar to resu l ts rep orted by [Marietta et al.l 
(|2000l ) and IPakmor et all (|200l for the mam-sequence 
case, but with a more compact companion and smaller 
binary separation. As a result of the asymmetric interac- 
tion, a small kick velocity, Wkick 85 km s^^, is imparted 
to the helium star. 

3.2. Parameter survey 

^ Movies are available at 
|http : / / sipapu . astro . illinois . edu/f oswiki/bin/view/Main/BinarySupernQvae| 



We conducted a parameter survey to explore the de- 
pendence of the numerical r esults on the bmary progen- 
itor's properties. Note that iWang et al.l ()2009bD deter- 
mined the orbital period, and thus the orbital separa- 
tion, at the onset of the supernova explosion, for each 
choice of helium star and white dwarf. In our study we 
do not follow the binary evolution up to the explosion, so 
the Roche- lobe radii of Wang et al.'s models are actually 
larger than the radii of our helium-star models. There- 
fore, in addition to varying the mass of the helium star, 
we also vary the binary separation in order to determine 
the effect of this parameter. 

Figure [3] shows the typical time evolution of the 
amount of mass removed from the helium star (defined as 
total unbound helium mass) by the SN ejeeta. The first 
peak occurs at the initial impact (~ 30 s) , but then some 
gas becomes bound again after the ejecta pass through 
the helium star. At ^ 80—100 seconds, a second peak de- 
velops which is associated with the effect of the reverse 
shock. Lastly, a third peak at ^ 500 seconds occurs 
when the compressed helium star relaxes and starts to 
oscillate. After 1,500 seconds, the amount of unbound 
mass reaches an approximately steady value. The final 
unbound mass is thus calculated for a given initial bi- 
nary separation and helium-star model by averaging the 
unbound mass values computed for several time steps at 
this late stage. 

Figures [4] and [5] show the final unbound mass and 
helium-star kick velocity (defined as the center-of-mass 
velocity of the final bound helium) as functions of the 
binary separation in our simulations. We include 14% 
error bars in Figure S] and Figure [5] based on the results 
of our convergence test (ii l3.3[) . We find that the unbound 
mass can be fit by the relation 



(1) 



where a is the orbital separation, mub is the power- 
law index, and the constant Cub depends only on the 
helium-star model (see Table[T|). For comparison, we also 
plot the powe r -law r elation with index —3.49 found by 
iPakrnor et al.l (|2008| ) and the data from iMarietta et al.l 
(|2000f ) for main-sequence companions (consistent with an 
index of —3.14). The power-law indices for our helium- 
star companions vary in a small range and bracket their 
results, suggesting that the index may be insensitive to 
the evolutionary state of the companion. The normaliza- 
tion of the above relation does appear to be sensitive to 
the nature of the companion star. 

The kick velocity of the helium star is calculated by dif- 
ferencing the center-of-mass positions in the z-direction 
at different timesteps. Alternatively, it can be deter- 
mined by dividing the total bound helium star momen- 
tum by the total bound helium star mass. Both methods 
yield the same kick velocity. To smooth out short-term 
fluctuations, we average the center-of-mass positions for 
every ten steps, then determine the kick velocity by nu- 
merical differentiation of these averaged values. The he- 
lium star is initially accelerated by the SN ejecta with an 
acceleration ~ GM^d/a^ for a time ~ a/wsn, and then 
it is kicked by the SN ejecta during the initial impact. 
During the initial impact, the kick velocity varies dra- 
matically. However, after ~ 1,000 seconds, these varia- 
tions decay. Thus, we deflne the kick velocity by using 
the difference between the final averaged velocity and the 
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Fig. 2. — The density distribution for model He-WDc (sec Tabic [T} with initial binary separation of 4 X 10^" cm. 
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Fig. 3. — Unbound helium mass versus simulation time using 
different levels of refinement for the He-WDc model (see Table [J| 



with initial separation 3 X 10 



(Note: the 10/12 run was very 



expensive and was carried out only up to 292 s.) 

maximum velocity just before the impaet of the SN ejecta 
(the exact time range to average the final velocity varies 



from run to run and covers the entire period of "smooth" 
variation in the kick velocity) . Figure [5] shows the kick 
velocity for the different helium-star models and initial 
binary separations. For initial binary separations larger 
than 4 i?Hc the kick velocity could not be adequately de- 
termined, because the perturbed velocity from the SN is 
larger tha n the kick veloc i ty at larger separations . As ob - 
tained by iPakmor et all (I2008D and iMeng et al.l ()2007D . 
a power-law relation is also found in our simulation and 
can be fitted by the relation 



kickfl 



(2) 



where i^kick is the kick velocity, m^ick is the power-law in- 
dex, and the constant Ckick depends only on the helium- 
star model (see Table [1]). However, unlike the situa- 
tion for the final unbound mass, the slope is very differ- 
ent from that found in the case of main-sequence com- 
panions (TOkick = — f .45 in IPakmor et al.l ()2008f ). and 
-1.26 in [Marietta et al.l (|2000D rfsee §4.1). 

3.3. Convergence test 

In order to determine the robustness of our numeri- 
cal results, we performed a convergence test for model 
He-WDc. We carried out simulations with the same ini- 
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Fig. 4. — Final unbound stellar mass versus binary separation 
for different helium-star models. The separation is expressed in 
units of the helium-star radius. The lines show power-law rela- 
tions from numerical simulations for different helium-star models 
(Table [ij- Error bars are based on the 14% error determined in 
our convergence test. T he da sh-dot line is the fitted curve pro- 
vided by Pakmor ct al] 1)20081 ) and the circles are the data from 
[Marietta et al.i (|200Q). 
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Fig. 5. — Similar to Figure |4] but for kick velocity versus binary 
separation for different helium-star models. Error bars are based 
on the 4% error determined in our convergence test. 

tial binary separation (a = 3 x 10^" cm) but several dif- 
ferent maximum AMR levels, computing the amount of 
unbound mass in each run as a function of time. The 
unbound mass is found by calculating the difference be- 
tween the initial helium-star mass and the measured total 
bound helium mass at each step. The total bound he- 
lium mass is the sum of the helium masses in all zones for 
which the total energy is negative. The results, which are 
plotted in Figure [3l show that the final unbound mass 
lies in the range between the 8/10 and 10/12 runs. The 
unbound mass is sensitive to the level of turbulence that 
occurs near the surface of the helium star, particularly for 
resolutions higher than 8/10 (see Fig. [S]), so rather than 
observing normal convergence behavior, we find that the 
unbound mass for runs with different resolutions fiuctu- 



ates within a small range of values. Similar behavior is 
also found with other initial binary separations. Thus, 
the difference in unbound helium mass between the 8/10 
and 9/11 calculations (~ 0.01 Mq) is used to estimate a 
relative error. Because the unbound helium mass reaches 
an approximately steady value at late times, this differ- 
ence suggests an estimate of about ±14% for the relative 
error in our runs, since it is not clear whether the 7/9 
run is always the lower limit of unbound mass for dif- 
ferent models or different separations. Similar analysis 
gives a ±4% relative error in kick velocity. To allow for 
a feasible parameter study, a resolution of 7/9 levels was 
chosen for the whole suite of twenty runs. 

4. DISCUSSION 
4.1. Stripped and ablated mass 

The impact of SN la ejecta on a helium-star compan- 
ion is not as dramatic as for a main-se quence or red- 
giant companion like those considered by [Marietta et al.l 
(|2000f ) because helium-star companions are more com- 
pact. This can be seen in Figure [71 which shows the 
gas pressure profile along the z-axis at several different 
times during the evolution of model He-WDc with an 
initial binary separation of 4 x 10^'^ cm. While the pres- 
sure immediately behind the ejecta shock front is initially 
~ 10^^ dyne cm~^. by the time the front reaches the he- 
lium star it has dropped to ^ 10^'^ dyne cm~^. The 
shock is considerably weakened by the time it reaches 
the deep interior of the star. As a result, the amount of 
unbound mass is much lower than in the case for main- 
sequence companions for a given ratio of separation to 
radius, a/R ~ 3, even though the helium-star channel is 
characterized by smaller binary separations. 

If we assume that the WD accretion and subsequent 
explosion in our models are driven by Roche lobe over- 
fiow (RLOF), the physically appropriate initial binary 
separations for models a to d should be 3.11, 3.00, 2.84, 
and 2.72 times i?Hc- The total unbound helium masses 
for these models at these radii (Figure |4]) arc ~ 1 — 3% of 
the initial helium-star masses, which is consistent with 
the sugg e sted u pper limit of ^ 0.01 Mq determined by 
[Leonard! ()2007[) from observations. 

The unbinding of mass from the helium star results ei- 
ther from ablation (heating) or from stripping (momen- 
tum transfer) by the ejecta. To estimate the relative 
contribution of these processes, we examined the mixing 
of nickel in the helium-star material. Because stripping 
involves the physical displacement of gas from the he- 
lium star by the ejecta, this process is associated with 
the contact discontinuity between the two fluids. This 
discontinuity is unstable, and the resulting instabilities 
lead to mixing, which while numerical in origin, nev- 
ertheless signals contact between the nickel-rich ejecta 
and helium-rich stellar envelope. In contrast, ablation 
proceeds through the shock heating of envelope material 
ahead of the contact discontinuity, and because the shock 
is more stable, much less mixing is expected where ab- 
lation is dominant. The averaged ablated and stripped 
mass after the initial impact for the different cases stud- 
ied is shown in Figure [51 

We find that the ablated mass is comparable to the 
stripped mass for smaller binary separations, and the 
amount ablated is sensitive to the binary separation dur- 
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Fig. 6. — The density distribution for model He-WDc with initial binary separation 3 X 10^*^ cm. Each frame shows the evolution time 
at around one minute for different AMR levels. The color scale indicates the logarithm of the gas density in g cm~'^. 




Fig. 7. — Pressure profile along the 2-axis at different times for 
model He-WDc with an initial binary separation of 4 X 10^*^ cm 
and 9/11 levels of refinement. The solid line is the original pressure 
profile before the SN explosion. 

ing the initial phase (less than 100 seconds for the case of 
model He-WDc). The final results suggest that in most 
cases the amount of stripped mass can be as much as an 
order of magnitude greater than the ablated mass, de- 
pending on the binary separation. The ratio cannot be 
determined accurately because we did not trace individ- 
ual fluid elements in these simulations. 

A simple an alytical method has been used by 
IMeng at al.l (I2007i ) to estimate the amount of unbound 
matter based on conservation of momentum, ignoring 
the shock dynamics and the effect of ablation. This 
method yields a shallower power-law slope (-1.9) for un- 
bound mass versus separation than the range that we 
observe (-3.1 to -4). However, the power-law slopes for 
un bound mass a nd ki ck velocity versus separation found 
by IMeng et al.l ()2007| ) are similar to the slopes inferred 
from the hydrogen-rich stars studied by [Marietta et al 



Although the method adopted bv IMeng et al 



2003) is oversimplified, their result suggests that abla- 
tion in main-sequence binary companions may not be 
as important as for helium-star binary companions. Fig- 
ure [8] shows that in our simulations the ablated mass cor- 
responds to ~ 1 — 20% of the total unbound mass, with 
a greater fraction ablated for smaller orbital separations. 
Thus, the total momentum imparted to the helium star 
results from contributions by both direct impact of the 
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Fig. 8. — Comparison of averaged ablated and stripped helium 
mass after the initial supernova impact versus the binary sepa- 
ration for model He-WDc (upper panel). Error bars represent the 
standard deviation during the average time range. The lower panel 
shows the ratio of ablated to stripped helium mass. 

SN ejecta and shock heating, whereas for hydrogen-rich 
companions the former effect is a more important con- 
tributor than the latter. The difference in the relative 
contribution of these two effects for different types of 
companion star may explain the difference in the power- 
law indices for kick velocity. 

4.2. Nickel contamination 

The companion star can be contaminated by the mix- 
ing of supernova ejecta with the helium-rich material 
in its envelope, perhaps resulting in a detectable en- 
hanced iron abundance after the nickel radioactively de- 
cays. Since ^^Ni is used as a tracer in our model, we 
have determined the amount of nickel bound to the rem- 
nant helium star. In general, there is a tendency for a 
higher level of contamination for more massive helium- 
star companions or for larger orbital separations. In the 
former case, a massive companion presents a larger cross 
section for capture of supernova ejecta, whereas in the 
latter case, the ram pressure from the ejecta is lower 
for larger binary separations, resulting in a reduction of 
the amount of contaminated matter that is stripped off. 
However, the level of contamination in our simulations 
is highly dependent on the development of small-scale 
fluid instabilities at multifluid interfaces. Although we 
find little evidence for a simple relationship between the 
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amount of contamination and the nature of the hehum- 
star companion or the orbital separation, we can estimate 
an upper limit on the nickel captured by the helium-rich 
companion by identifying the ejecta which cannot escape 
the gravitational potential. This limit is ^ 5 x 10~^ Mq, 
somewhat smaller than the estimate by Marietta et al. 
(~ 1.5 X 10~^ Mq). If the mixing of nickel is only re- 
stricted to the envelope, the nickel to helium ratio can be 
estimated by using the upper limit of nickel contamina- 
tion and the final envelope mass. We define the core mass 
for the initial helium-star models by finding the radius at 
which the second derivative of the density with respect 
to radius is a maximum; the envelope mass is then the 
total mass less the core mass. Using this definition, we 
find that the unbound helium-star material is taken en- 
tirely from the envelope. The final envelope mass can 
thus be obtained from the difference between the initial 
envelope mass and the final unbound mass in Figure |4] at 
the separation corresponding to RLOF. We find that the 
ratio of the upper limit on bound nickel mass to the final 
envelope helium mass is ~ 9 — 50 x 10~^ for the different 
helium-star models. This value is substantially higher 
than the solar ratio of iron abundan ce to that of hydro- 
gen pl us helium (5.1 x 10"'') found by I Anders fc Grevessd 
(|1989D . suggesting that the abundance of nickel/iron in 
the remnant helium-star atmosphere should be enhanced 
relative to normal Population I stars if surface convection 
is not an important factor. 

4.3. Detecting the remnant companion star 

After the supernova, the companion star moves with 
a velocity v ~ Vo,-i, + VMck- The orbital speeds corre- 
sponding to the range of binary separations in our simu- 
lations are Voih ^ 350 — 800 km s^^, while the kick speed 
'^kick ~ 30 — 100 km s~^. Thus, the net velocity is pri- 
marily determined by the orbital motion for most binary 
separations. Figure |9] shows the ratio of kick speed to 
the total speed, assuming that the kick velocity is per- 
pendicular to the orbital velocity. If the mass transfer in 
the binary system is via RLOF, then the kick contributes 
little to the total velocity. For very small orbital separa- 
tions, the kick contributes at most ~ 10 — 20% of the 
total. 

Based on these observations, we may expect that he- 
lium stars with high space motion found to be associ- 
ated with Type la supernova remnants would be evidence 
for SNe la proceeding through the helium-star single- 
degenerate channel. If such objects are also found to have 
enhanced iron abundances relative to helium, this evi- 
dence would be considerably strengthened. In such cases 
the relative velocity between the helium star and the cen- 
troid of the supernova remnant would place a stronger 
constraint on the initial binary separation than on the 
initial helium star's mass or the white dwarf mass. 

5. CONCLUSIONS 

We have investigated the impact of SN la ejecta on 
a companion star in the single-degenerate helium-star 
channel for the short-delay-time population of SN la 
via two-dimensional hydrodynamical simulations using 
a range of helium-star models and binary orbital sep- 
arations. Although the general behavior of the tem- 



po ral evolution is simila r to that in previous stu dies 
by [Marietta et"all (|2000f ) and iPakmor et all (I2008D for 
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Fig. 9. — The ratio of kick speed to total speed for different 
helium-star models as a function of initial separation. Lines: from 
the fitted kick velocity. Data: from the calculated kick velocity. 
Error bars are based on the 4% error determined in our convergence 
test. We assume the kick velocity to be perpendicular to the orbital 
velocity. 

main-sequence and red-giant companions, the amount of 
matter unbound from the helium stars is less than for 
hydrogen-rich companions. Due to the shorter orbital 
periods of helium-star progenitor systems, the space mo- 
tion of the companion star after the explosion is found 
to be higher. We find a power-law relation between the 
unbound mass and initial binary separation that is con- 
sistent with the previous studies, suggesting that the 
power-law behavior is not strongly sensitive to the na- 
ture of the companion. The kick velocity can also fitted 
by a power-law and we conclude that the power-law in- 
dex may refiect the relative importance of the effect of 
ablation. An upper limit on the amount of nickel cap- 
tured by the helium star is found to be '--^ 5 x 10"'' Mq. 
The ratio of nickel to helium abundance may be useful as 
a diagnostic of such events in future observational studies 
of SN la stellar remnants. 

Future work in this area will include relaxing the as- 
sumption of axisymmetry to model the mixing of ejecta 
with the helium star by including the binary orbital mo- 
tion in three spatial dimensions. Additionally, including 
radiation transfer within these simulations will allow us 
to determine how much of the helium is ionized, allowing 
us to make direct contact with spectroscopic constraints 
on the presence of helium in these systems. 

We thank the anonymous referee for his/her valuable 
comments and suggestions. The simulations presented 
here were carried out using the NSF Teragrid's Ranger 
system at the Texas Advanced Computing Center un- 
der allocation TG-AST040034N. FLASH was developed 
largely by the DOE-supported ASC/ Alliances Center for 
Astrophysical Thermonuclear Flashes at the University 
of Chicago. This work was supported, in part, by NSF 
AST-0703950 to Northwestern University. 
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